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ABSTRACT 

We present TA-DA, a new software aimed at greatly simplify and improve the analysis of stellar 
photometric data in comparison with theoretical models, and allow the derivation of stellar parameters 
from multi-band photometry. Rs flexibility allows one to address a number of such problems: from the 
interpolation of stellar models, or sets of stellar physical parameters in general, to the computation 
of synthetic photometry in arbitrary filters or units; from the analysis of observed color-magnitude 
diagrams, to a Bayesian derivation of stellar parameters (and extinction) based on multi-band data. 
TA-DA is available as a pre-compiled IDL widget-based application; its graphical user interface makes 
it considerably user-friendly. In this paper we describe the software and its functionalities. 
Subject headings: Data Analysis and Techniques, Astrophysical Data, Stars, Star Clusters and Asso- 
ciations 



1. INTRODUCTION 

A common problem in stellar astrophysics is the rela- 
tion between the observational quantities (magnitudes, 
colors, fluxes) and the physical parameters of the stars 
(effective temperature, luminosity, mass, age). Often, 
dust extinction challenges this task even further, espe- 
cially when differential reddening affects the individual 
sources in a region. Whereas it has always been common 
practice to analyze photometry data in 2-dimensional 
space, like a color-magnitude diagram (CMD) or a 2- 
color diagram (2CD), often this approach is insuffi- 
cient to fully characterize the individual stars. In fact, 
there are well-known degeneracies among the stellar 
parameters for stars of given color: a few examples 
are the effective temperature (Tes) -extinction degener- 
acy for low- and intermediate-mass stars and optica l 
magnitudes (e.g.. iHillenbrandl [l997l: lBailer-Jonesll2011[ ). 
or the age-metallicity degeneracy for RGB s tars (e.g., 
iVandenBerg et al.|[200l iCrnoievic et al.ll2010[ ). 

Multi-band photometry often allows to improve upon 
these limitations, both in breaking the degeneracies 
and reducing the overall uncertainty. Today, the avail- 
ability of large field s urveys such as the S loan Digi- 
tal Sky Survey (SDSS, Abazaiian et al. '2009) , the Two - 
Micron All-Sky Survey (2MASS, Skrutskic e taIll2006D 
and the Wide-fiel d Infrared Survey Explorer (WISE, 
iWright et al.ll2010l ). the diffusion of pubhc data archives 
and the rise of Virtual Observatory (VO) infrastructures, 
enable easy access of an overwhelming quantity of pho- 
tometric data for the community. 

Even with a great variety of data in one's hand, the 
characterization of the properties of stars may hide fur- 
ther obstacles, some of which merely practical. Various 
theoretical models predicting the stellar parameters are 
known to differ from one another, requiring one to as- 
sess, by comparison with the data, the level of trustwor- 
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thiness of the models. Often, stellar evolutionary models 
already provide the predicted magnitudes and colors; this 
conversion, however, may not be well calibrated, or not 
available in the specific photometric bands that one re- 
quires. In this case, common solutions are to transform 
the available photometry from one to another system - 
with the consequent uncertainties - or to rely on syn- 
thetic photometry to recompute the predicted fluxes. 

There are several available tools to simplify some of 
these steps. A number of stellar evolutionary mod- 
els provide online interfaces for interpolation of their 
grids and compu tation of magnitudes and colors (e.g, the 
ISiess et al.l[2000l pre-main sequen ce models or the Padov a 
models from lMarigo et al.ll2Ci08l and lGirardi et al.ll2010[) . 
Similarly, some grids of synthetic spectra already pro- 
vide, besides the computed spectra, also tools to compute 
the integrated colors and magnitudes in several photo- 
metric syst ems (e.g., the PHO ENIX models of the Lyon 
group, see lAllard et all 1201111 . For what concern the 
synthetic photometry in general, a very popular code is 
SYNPHDT, an IRAF package part of the Space Telescope 
Science Data Analysis System (STSDAS) which per- 
forms several different types of such calculations. As for 
the estimate of stellar parameters and extinctions from 
multi-band photomet ry, a powerful code is CHORIZOS 
(jMai'z- Apellanid [2004D . This software allows one to es- 
timate Teff, logg, the extinction Ay as well as the red- 
dening parameter Ry from multi-band colors, using a 
Bayesian approach. CHORIZOS also performs automatic 
synthetic photometry, with the possibility of adding cus- 
tom filter bands. Still, CHORIZOS is designed to fit stellar 
colors, not magnitudes, neglecting the additional infor- 
mation in the data provided by the luminosities; more- 
over, stellar evolutionary models are not directly inte- 
grated in the code. All these tools, in fact, have the com- 
mon shortcoming that either they are designed to address 
only a particular problem within the general analysis, or 
that they are constrained to perform as black-boxes. 
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Fig. 1. — Some examples of the graphical user interface of TA-DA: top left: the interpolation of evolutionary models or the definition of 
the model stellar parameters; top right: the interface for the synthetic photometry; bottom right: the visual inspection of the computed 
synthetic photometry, the attachment of observed photometry and the selection of sources within the observed catalog; bottom left: the 
interface for the stellar parameter fitter. 



Following these premises, it appears that a new, in- 
tegrated, flexible code for photometric data analysis in 
relation to stellar parameters is needed. With this mo- 
tivation, we have developed the Tool for Astrophysical 
Data Analysis (TA-DA). In this paper we present the 
code, its general functionalities, and show some example 
of use. 

2. REQUIREMENTS FOR TA-DA 

TA-DA is conceived as a tool for quick-look, compari- 
son, and parameter fitting of multi-band photometry in 
relation to stellar models. In the development of the 
code, we focused on two general concepts: 

1. Generality and Versatility: TA-DA should be able 
to address a range of different problems concern- 
ing the relation between photometry and stellar 
parameters. 

2. Usability: TA-DA must be user-friendly, ensuring 
the users to exploit its functionalities (or test if it 



is useful for their scientific aims) with little or no 
extra time dedicated to familiarize with the code. 

Specifically, we included the following functionalities: 

• Integration with stellar evolutionary models: TA- 
DA is natively able to load and interpolate stel- 
lar interior models, allowing one to fit model- 
dependent quantities such as stellar masses and 
ages. The consideration of interior models, never- 
theless, is not mandatory, as ultimately one could 
solely consider, e.g., Tcff as the only stellar param- 
eter. 

• Automatic and universal synthetic photometric en- 
gine: TA-DA is able to perform synthetic photom- 
etry, converting stellar parameters into magnitudes 
or fluxes in a number of units, allowing also redden- 
ing due to dust extinction to be directly applied to 
the synthetic spectra. The synthetic photometry 
can be performed on different grids of atmosphere 
models, and on arbitrary photometric systems. 
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• Data visual quick-look in comparison with models 
or theoretical stellar parameters: the user can load 
photometric tables, and plot them together with 
the pre-computed sets of stellar parameters, con- 
structing, e.g., CMDs or 2-CDs, and varying the 
amount of assumed dust extinction. This is useful 
for obtaining a first visual check on the behavior of 
the data, and refine the selection of sources to be 
further studied. 

• Automatic parameter fitter: TA-DA allows one to 
automatically estimate the stellar parameters of 
the individual sources from their multi-band pho- 
tometry; depending on the data at disposal and 
the number of free parameters, this can be either 
an exact solution or a probabilistic one. The num- 
ber of free parameters, as well as the photometric 
bands to use, can be fully customized. 

3. HOW TA-DA WORKS. 

TA-DA has been developed using the ITTVIS Inter- 
active Data Language (IDL), and presents a graphical 
user interface (GUI) based on widgets; it is distributed 
as a pre-compiled, cross-platform code. The distribution 
package already includes a number of stellar evolution- 
ary and atmosphere models, as well as a full repository 
of photometric filter throughput profiles. We refer the 
readers to the Appendix [K\ for additional details. 

The TA-DA main GUI is divided into 4 panels, re- 
spectively for: a) the specification of the model stellar 
parameters and the interpolation of evolutionary mod- 
els; b) the synthetic photometry; c) the visual display of 
the synthetic photometry, the handling of the observed 
photometry table, selection and plotting of the latter; d) 
the stellar parameter fitter. We will now describe some 
of these functionalities. 

3.1. Stellar parameters and evolutionary models 

TA-DA includes several families of evolutionary 
models; for pre-main sequence (PMS) popu lation s 
we have included models from Baraff e et al 
D'Antona &: Mazzitellil (I1998D : [PaUa fc Stable: 



1998); 
1991 : 



Siessetal.1 (|2000[) : iTognelh et all (|2011[ ). for evolved 
population we inc luded th e Padova ni odels from 
iMarigo et all (|2008D with the iGirardi et all (2010) cor- 
rection. If multiple values of metallicity are available for 
a given set of interior models, a value should be speci- 
fied. When instead a grid of models includes additional 
parameters, such as in the case of the PISA models from 
iTognelli et al.l ()2011[ ). computed for several values of mix- 
ing length parameter a, helium enrichment, and deu- 
terium initial abundances, these are treated as different 
families of models. Additional interior models can be 
added by a user provided they are stored in the proper 
format. 

TA-DA allows one to consider the entire grid, for a 
selected metallicity, or to restrict it to a particular arbi- 
trary isochrone or mass track. Furthermore, it is possible 
to manually input user defined sets of stellar parameters, 
attaching a table reporting an arbitrary set of 2 quanti- 
ties chosen from the following 5: mass, log age, log Toff, 
logl/boi, stellar radius. It should be noted that any com- 
bination of 2 of these parameters automatically define the 
remaining ones, since, e.g, the stellar radius R relates 



to Tcff and logL through L = AnR^T^g and the evo- 
lutionary models natively map the Hertzsprung-Russel 
Diagram (HRD) into stellar masses and ages. Thus TA- 
DA will consider the user-supplied parameters and derive 
the remaining ones, including the surface gravity logg, 
needed for the synthetic photometry (see below). 

Alternatively, TA-DA also allows one to consider only 
parameters in the HRD (Tcs, R, logg for a specified 
metallicity, if available), the minimum ones necessary 
to both derive the synthetic photometry and allow the 
fitting of stellar parameters, without relying on model- 
dependent physical quantities such as stellar masses and 
ages. 



3.2. The synthetic photometry 

TA-DA is able to perform automatic synthetic pho- 
tometry, i.e., to convert stellar parameters into observed 
fluxes or magnitudes in a given photometric filter, rely- 
ing on grids of synthetic spectra. To this purpose, one 
must select a grid of atmosphere models and, if avail- 
able, the metallicity [M/H]. Dust extinction can be also 
added (by declaring one or several values of Ay), and 
directly applied to the synthetic spectra before the com- 
putation of the synthetic phot ometry. We in c luded the 
Galactic reddenin g laws from iCardelh et all (|1989f ). as 
well as those from iGordon et al.l (|2003f ) fo r the LMC av- 
erage, LMC2 supershell and SMC bar from lGordon et all 
l|20M ). In first case, besides Ay one (or more) values 
of the reddening parameter Ry must be specified. The 
distance modulus in magnitudes must also be provided. 
Finally, one specifies the photometric system and bands 
in which the synthetic photometry must be performed 
(see the Appendix [Al for details). 

Synthetic photometry is performed in the standard 
way, by integrating the synthetic spectra within the fil- 
ter bandwidths and normalizing onto a spectrum of Vega; 
the observed magnitude in a photometric band S\ of a 
star with a spectral energy distribution F\ , stellar radius 
R, and true distance modulus /j, = (m — M)o is given by 



Ms 



-2.5 log 



R 
lOpc 



XfxSxdX 



(1) 



where is a reference spectrum that gives a known 
apparent magnitude ZPs^ ; in the Vegamag standard, 
which uses the fiux of a Lyr as reference, — -F'a,vega 
and the zero-points ZPs^ should be zero by definition, 
although this is not always the case for some old photo- 
metric systems. Finally, A\ is the reddening law, scaled 
to the selected values of Ay. Equation ([T]) can be rewrit- 
ten then as: 



Ms, =-5 log 



(^)-51og(fe)+^^ + 5(^A,5.) 
43.2337- 5 log ) +^ + B(Fa, 5a) (2) 



4 



Da Rio et al. 2012 



where 

B{Fx,Sx) ^-2.5log(^ (3) 

The latter term can be directly calculated having the syn- 
thetic spectrum for the Tea and log g corresponding to a 
given point of the input model parameters, a calibrated 
Vega spectrum and the band profile. The Vega normal- 
ization is performed using a recent reference spectrum of 
Vega (Bohlin 2007). We have co nsidered the zeropoints 
ZPs^ from.Maiz-Apellaniz (' 2007D . collected from several 
sources for the code CHORIZDS. 

Although the default units for the computed syn- 
thetic photometry are Vegamag magnitudes, TA-DA al- 
lows to switch also to ABmag, STmag, Jansky or erg 
s~^cm~^A~^. This is performed by computing, for ev- 
ery filter, the flux corresponding to m = mag by direct 
integration onto the spectrum of Vega. Then, since AB- 
mag and STmag are related to constant flux densities 
(respectively per unit of spectral frequency and wave- 
length), these "zero-point" fluxes automatically provide 
the offsets between the 3 magnitude standards. 

TA-DA allows one to save the results of synthetic pho- 
tometry to file, as well as to use them for the subse- 
quent analysis (visual inspection and parameter fitting). 
One can also save a table with most of the relevant 
characteristics of the selected filters. These include, be- 
sides the aforementioned zeropoints, also effective wave- 
length, full width half maximum (FWHM) and equiv- 
alent width (EW). This information can be particularly 
useful when one adds new (non-conventional and not pre- 
characterized) filters to the tool, as TA-DA automatically 
derives these quantities on the fiy 

TA-DA plots the results of synthetic photometry on 
screen on a dedicated window. Here one chooses the 
quantities to be plot, either a magnitude-magnitude 
plane, a CMD, or a 2CD. The units and extinction to 
be considered can also be edited in real time. An inter- 
esting feature that is helpful for a fast quick look at the 
result of synthetic photometry is that one can see the 
parameters associated to individual points of the model 
grid by clicking directly on the plotting window. All the 
physical parameters associated to that point, as well as 
the displayed observational quantity (e.g., color and mag- 
nitude) are then reported in a separate popup window. 

3.3. Observed photometry vs. models 

TA-DA enables one to load a table reporting observed 
photometry, both for quick look in comparison with the 
models and for the fit of the stellar parameters. Cur- 
rently the software supports both ASCII formats and 
Virtual Observatory table. A new GUI window is cre- 
ated, reporting attached table, and enabling one to 
browse through the data, sort by columns, and select 
individual sources. In this panel, one must specify what 
columns correspond to what filter, as well as the associ- 
ated uncertainties (if available); the options are chosen 
among the photometric bands previously selected for the 
synthetic photometry. The selection of individual sources 
can be performed both manually (clicking on table en- 
tries), or automatically, through logical conditions not 



restricted to the photometric data. An example could be: 
select all the source with V < 20, with the photometric 
errors al < 0.1 and where the last column of the table 
does not report the string 'Class F. The selected sources 
are also automatically shown in the plotting window of 
TA-DA, together with the models previously computed. 
Similarly as for the models, one can click directly on in- 
dividual stars on the plot, and see on screen the location 
of that source in the catalog, as well as the value for that 
star of all the columns in the table. 

3.4. The parameter fitter 

With all the necessary preparation completed (a cat- 
alog of observed magnitudes or fluxes and an a grid of 
models converted into the same observables through syn- 
thetic photometry) , TA-DA is able to perform a series of 
fitting techniques to estimate the parameters of the ob- 
served sources. Specifically, it allows one to choose which 
bands to consider, and if magnitudes or colors are to be 
used. The latter options is useful when one aims to de- 
rive the stellar parameters in a luminosity (or radius) 
independent fashion, in a 2-color or multi-color space. 
Clearly, choosing colors instead of luminosities decreases 
the dimension of the observational space by one unit. 

As for the theoretical space - the result of synthetic 
photometry - in general one could leave from 1 to 3 free 
parameters to be evaluated by the fitter. One or two pa- 
rameters correspond to the photosphere. This depends 
on the model previously assumed for synthetic photom- 
etry: for instance, if one previously considered 1 single 
isochrone, or mass track, or arbitrary curve in the (Toff- 
log g) plane, etc., this will lead to a 1 free dimension 
in the parameter space. Conversely, considering an en- 
tire grid of interior models (masses and ages), or a 2- 
dimensional grid of Toff and stellar radii brings the free 
photospheric parameters to be determined to two. The 
additional free parameter is the extinction Ay for indi- 
vidual stars. TA-DA allows to leave Ay unconstrained 
within a pre-defined range, as well as to constrain it to 
a given value for all stars. If synthetic photometry was 
previously carried out for multiple values of the redden- 
ing parameter Ry, which is possible if the lCardelli et ahl 
p989i) reddening law is assumed, the value of Ry to be 
used by the fitter must be specified. We did not al- 
low to leave both Ay and Ry free, as for example in 
the tool CHORIZDS by Mafz-ApcUaniz (2004). This is 
because, in general, most photometric datasets are not 
accurate enough to allow a precise disentanglement of 
the two; moreover inaccuracies in the synthetic spectra 
would likely induce large deviations in the derived values 
of Ry. 

In order to make the fit possible, the dimension of the 
observational space m (number of colors or magnitudes) 
must be greater or equal than the number of free param- 
eters n. Although from the algorithmic point of view 
there is no separation between the cases m > n and 
n = m, there are some conceptual differences. The first 
case is analogous to a standard SED fitting, in which the 
model space is explored (and if needed interpolated to 
a finer sampling) until the minimum is determined 
(see Equation m below) ; if x^j„ — 1, the best solution 
is also a good fit, otherwise, the data may be incom- 
patible with the models, or the associated uncertainties 
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underestimated. The second case, on the other hand, 
corresponds to an interpolation problem, for which the 
solution is expected to have = 0; if > Oj the lo- 
cation of the data point in the observational spaces lies 
outside the area (or volume) covered by the models. 

In detail, the fitter of TA-DA works as follows: let us 
consider P the rt-dimensional space of the free param- 
eters (pi,...,p„); the synthetic photometry associates to 
each point in P a point (si, Sm) in the m-dimensional 
space of the observational quantities, S. Then, for ev- 
ery observed star, the observed fluxes, colors or mag- 
nitudes (ri,...,rm) with the associated uncertainties 
((Tri, fr^) are considered and the is simply: 

X^bi,...,p„)=f;ill^. (4) 

First TA-DA identifies the set of (pi,...,p„) which mini- 
mize x^, then P is locally oversampled to a denser grid 
P2, and S is interpolated accordingly to ^2. The x^ is 
thus recomputed, and a new set of parameters estimated. 
The process is iterated until a convergence is found, ei- 
ther ~ for n = m or to a positive value otherwise. The 
local oversampling of P is made easier by the fact that 
the grid of parameters is always rectilinear in (pi,...,p„). 
This could be generally, e.g., a rectilinear grid of masses, 
ages and Ay or any other choice for the original model 
stellar parameters, for instance a rectilinear grid of Teff, 
stellar radii i?, and Ay- Since, as mentioned earlier, the 
stellar parameters mass, age. Teg, logL, and R are inter- 
connected, for each solution all these quantities will be 
provided, irrespective of the type of initial parameters 
provided to TA-DA. 

In the case n = m, as already mentioned, a x^ > 
indicates that the data point lies outside the volume 
covered by the models. In practice, the fact that fitter 
converges until the x^ does not improve by local over- 
sampling within P always leads to small positive values 
of x^- A source located outside the model grid, but by 
less than one unit of photometric error from the edge, 
will also lead to a x^ < 1- To precisely distinguish if 
the point is inside or outside the model grid, therefore, 
TA-DA also determines this condition geometrically, re- 
gardless the minimum x^ > 0. In any case, even for data 
points incompatible with the models, the closest solution 
is found, and the associated x^ tabulated in the result. 
The user must therefore check manually the results and 
decide if there are any solutions to be rejected, based on 
the specific scientific goals. 

The TA-DA fitter also allows one to estimate the un- 
certainty associated to each parameter of the best fit- 
solution for the individual stars. This is performed by 
TA-DA with a Monte Carlo simulation, in which the 
photometry is displaced according to the photometric er- 
rors. This approach may be computationally time con- 
suming, especially in the case of large grids of models 
and 3 free parameters, or large observational data sets. 
This method, however, remains the safest for an accu- 
rate estimate of the errors; this is because the relation 
between P and S provided by the synthetic photometry 
is often highly non-linear, and symmetric, well-behaved 
photometric errors in the observational space relate to 
highly correlated and skewed probability distributions in 



the model space. 

After the completion of the fitting procedure, TA-DA 
allows one to save the results to file. These include all 
the best-fit parameters, with the associated uncertain- 
ties if available, the minimum x^, and a flag parameter 
indicating if the solution is a good fit or if the observed 
photometry was incompatible with the models. 

4. IMPORTANT CAVEATS 

4.1. Accuracy of the synthetic spectra 

It should be beard in mind that the correctness of the 
fitted parameters is strongly related to the accuracy of 
the assumed models, e.g. both the evolutionary mod- 
els and the synthetic spectra. It is well known that 
different families of evolutionary models differ in their 
predicted masses and ages from eac h other, and this is 
particularly sever e for PMS models (iMavne et al.l [20071 : 
iHillenbrandl |2009() : there is also evidence of systematic 
offsets between the mod eled masses of low-mass star s and 
the measured once ('e.g.. lHillenbrand fc Whita[200l . On 
the other hand, current atmosphere models do not seem 
to fully reproduce the observed fluxes, especially for low- 
mass M-type stars and at young ages , both at optical 
wavelengths (e.g., iDa Rio et aT]|2009U 120101: iBell et al.l 
120121) and, although less prom inently, in the near in- 
frared (|Scandariato et al1l2012[ ). Because of these issues, 
a dataset compared with models in different colors can 
lead to systemati c differences in the derived parameters 
(e.g., iNavlor et al.ir2002i ). 

We auspicate that developments in the stellar atmo- 
sphere modeling will solve these problems; alternatively, 
it would be useful to produce empirically calibrated 
grids of spectra, constructed to match the observed col- 
ors. This has been attempted in the past for main- 
sequ ence dwarfs and gia nts (see, e.g., the BaSeL grid 
from lLeieune et aI]|1997D . but such a calibration is still 
missing for, e.g., PMS stars. 

In any case we warn the users of TA-DA to be aware 
of the accuracy of the synthetic spectra grids to be used 
in the parameter range relevant for their astrophysical 
analysis. We also stress that the plotting and quick look 
capability of TA-DA (e.g., looking at the same isochrones 
compared to the observed data in different color-color 
and color-magnitude planes) can be very helpful to pin- 
point possible problems with the fluxes predicted by the 
models, as well as possible errors in the photometric data 
calibration. 

4.2. N on- photo spheric fluxes 

As a probably unnecessary note, we highlight that TA- 
DA is meant for estimating stellar parameters in the case 
that the observed fiuxes depend only the photospheric 
emission and presence of dust extinction. When this is 
not true, such as when the fiuxes are contaminated by 
circumstellar emission or other fiux excesses, the results 
obtained could be wrong. In some cases, however, it 
is possible to use TA-DA for non standard applications, 
such as disentangling photospheric properties from flux 
excesses. As an example, having at disposal infrared 
multi-band photometry, one can isolate the excess due 
to circumstellar disks around young stars and derive T^e 
for individual stars. To do so, one could create a tuned 
grid of synthetic spectra in which instead of log g as sec- 
ond parameter (besides Teg), the amount of disk excess 
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Fig. 2. — PMS evolutionary tracks from ribgnelli et a l.l lj201U) converted into the ACS photometric system using TA-DA, together with 
the observed photometry of LH 95. Sources are color coded according to their masses (left panel) and ages (right panel), as derived by the 
TA-DA fitter. Open squares are sources that the software identifies to be outside the range covered by the models. 

appropriately the Ay range allowed during the fitting 

0.3^ — — — — ' — ■ — ■ — — — • — — — — — • — — — — — □ can also solve the problem. 
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Fig. 3. — Uncertainties in the LH 95 members masses and aged 
derived by TA-DA, as a function of these two quantities, originated 
from the individual photometric errors. 

is added. Then TA-DA can be run on a multi-color ob- 
servational space, and the fitter will provide the solution 
for both internal and external effects on the stellar flux. 

4.3. Degeneracies in the solutions 

In some cases one could encounter degeneracies be- 
tween the stellar parameters. For example, when dered- 
dening a CMD onto an isochrone, for some sources there 
could be multiple intersection between the de-reddening 
line and the model. Since the fitter of TA-DA simply 
minimizes the x^: in case of multiple solutions, only one 
is selected. However, the Monte-Carlo simulation used to 
estimate the uncertainties will likely converge on differ- 
ent solutions in different iterations; this leads to a larger 
associated uncertainty for the parameters, and this can 
be used as a diagnostic tool when analyzing the results. 
If multiple solutions are caused by different intersections 
of the de-reddening line with the models, constraining 



5. EXAMPLES OF APPLICATIONS OF TA-DA 

We present here some practical examples of applica- 
tions of TA-DA, for illustrative purposes. 

5.1. CMD interpolation 

Let's suppose we have a CMD and we aim at in- 
terpolating isochrones and tracks to estimate ages and 
mass es of the sou r ces. W e consider the photometry 
from iDa Rio et al.l (|2009aD . obtained with the Hubble 
Space Telescope (HST) Advanced Camera for Surveys 
(ACS) on the young cluster LH 95 in the LMC, using 
the fihers F555W an d F814W. We r un TA -DA, assum- 
ing the models from iTognelli et al.l (|2011[ ) with mixing 
length parameter a = 1.68, normal initial He and Deu- 
terium abundance, and the LMC metallicity Z = 0.008. 
We perform the synthetic photometry in our 2 filters 
with TA-DA cons idering the synthetic spectra BT-Settl 
from lAllardeFaLl pOll.) , for the same metallicity as for 
the interior models; we assume a true distance modulus 
fi = {m — M)o = 18.41, and two values of Ay'- and 0.5 
(the average extinction for L H 95), consider i ng the LMC 
average reddening law from IGordon et al.l (|2003D . We 
save the synthetic photometry to file for later use, and we 
load the table with the observed photometry to TA-DA. 
We run the fitter, constraining Ay = 0.5 and perform- 
ing the MC simulation for the uncertainties. Results are 
shown in Figure [2l here we show the models converted 
in our CMD plane, together with the sources color-coded 
according to their estimated parameters from TA-DA. It 
is evident how TA-DA correctly interpolates the models 
to the observations, deriving masses and ages. It also 
identifies the sources located outside the model grid, and 
assigns them the parameters of the closest model point. 
There is some scatter in the assigned ages for interme- 
diate mass stars: this is due to a residual degeneracy in 
the models in this parameter range, where a given point 
in the CMD can be associated both to a PMS and post- 
MS evolutionary stage. In Figure [3] we plot the errors 
in age and mass as a function of the derived parameters, 
due to the photometric errors, and computed by TA-DA 
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Fig. 4. — Same as Figure 3] for the second example of Section [S] Here we consider NIR photometry of the ONC, and leave as free 
parameters the photospheric locus on a 2 Myr isochrone, as well as the extinction Ay, constrained to be from to 10 mag. The data in 
the left and right panels are color-coded, respectively, for Tcff and Ay as derived by the TA-DA fitting routine. 

using a MC simulation. The age uncertainty increases j 

with ages, for the reason mentioned above. The increase 
in the mass uncertainty do the low-mass range, on the 
other hand, is due to larger photometric errors at faint 
luminosities. 

5.2. De-reddening to an isochrone 

In this second example, we consider a CMD and use 
TA-DA to deredden all the sources to a single given 
isochrone. We use the JHK photomet ry of the Orion 
Nebula Cluster (ONC) presented in IRobberto et al.l 
(jlOlO). We run TA-DA con sidering a single 2 Myr 
isochrone from the models of IBarafFe et al.l (jl998f) , we 
apply the distance modulus for the ONC of 8.085 mag, 
consider the BT-Settl models for solar metallicity, and 
obtain the synthetic photometry in the 2MASS photo- 
metric system (which is the one in which our data are 
expressed. We run the fitter leaving Ay as a free pa- 
rameter, in the range < Ay < 10 mag. The result is 
illustratively shown in Figure HI As for Figure [H open 
dots represent sources outside the considered range of pa- 
rameters; it can be noted that all the sources at red colors 
{J — H >2) appear as such. This is because of the upper 
limit imposed for Ay. We also highlight that for these 
sources, the parameters (indicated by the color-coding 
of the data points) represent the "closest" model points, 
along the direction defined by the individual photometric 
errors in the J, H 2-magnitude plane. This is the reason 
why, in the left panel of FigureS|and for {J — H > 2), the 
estimated Tcfi appears to decrease with increasing color 
(J-H) 

5.3. A multi-band example 

In the next illustrative example we show the capa- 
bilities of TA-DA in dealing with multi-band photom- 
etry. We start by simulating a population of 500 evolved 
sources, with random masses in the range 1 < M/Mq < 



1 

-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 

B-V 

Fig. 5. — Our simulated photometry of 500 stars with masses 
between 1 and 10 Mq and ages between 10 Myr and 1 Gyr, in 
the BV CMD. The upper panel shows the intrinsic magnitudes for 
Ay = 0, color-coded according to the stellar age. The lower panel 
shows the same population after applying a random differential 
extinction in the range < Ay < 2 and photometric errors. 

10 and random ages between 10 My r and 1 Gyr. We uti - 
lize TA-DA, assuming the models of lMarigo etall (|2008[ ) 
first to compute the photospheric parameters (Tos, ra- 
dius, etc.), then to derive by means of synthetic pho- 
tometry their magnitudes in the standard UBVRIJHK 
filters (Johnson, Cousins, and 2mass systems). An ex- 
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Fig. 6. — Comparison of the original parameters of the 500 simulated stars, and the parameters recovered by the fitter within TA-DA 
after adding a random amount of extinction and photometric errors. The three rows shows the result of the fit using only 3, 5 or all 8 
UBVRIJHK bands. 



2 




2 4^5 a 10 

Fig. 7. — The error in the estimated ages as a function of stellar 
mass. 

ample of the result, limited to the BV CMD, is shown in 
Figure [SJ upper panel. Next, we add a random amount 
of extinction to each of the 500 simulated stars, in the 
range < Ay < 2. We further displace the magnitudes 
by adding random, uncorrelated photometric errors to all 
the magnitudes, assuming an amplitude a — 0.05 mag. 
The result, again in the BV CMD is show in Figure [SI 
bottom panel; it is evident how at this point, a single 
CMD is insufScient to recover the properties of the indi- 
vidual sources, as they appear well mixed. 



We now use the fitter integrated in TA-DA to try 
to recover the original parameters of these 500 stars. 
We consider the e ntire grid of interior models from 
iMarigo et all (|2008[ ). convert them in absolute magni- 
tudes the UBVRIJHK bands, load the simulated ta- 
ble with the photometry of our simulated sample, and 
run the fitter to estimate the most likely parameters of 
the stars, leaving all mass, age, and Ay as free param- 
eters. We perform this in multiple configurations, vary- 
ing the number of photometric bands used to estimate 
the parameters. Some examples are shown in Figure 
[6l Here we compare the original mass, age and extinc- 
tion, with those automatically derived by the fitter, using 
only 3 bands (BVI), 5 bands [BVRIJ) or aU 8 bands 
{UBVRIJHK). It is evident that increasing the num- 
ber of bands leads to a better estimate of the parameters; 
with all 8 bands, TA-DA is able to fully disentangle the 
differential Ay from the stellar mass, thus recovering the 
properties of the sources. The residual scatter is due to 
the random photometric errors that we introduced. 

The central column of Figure [5] shows that, on the 
other hand, the estimate of the age remains highly un- 
certain. This is not a problem with the fitting procedure, 
since, as the extinction is correctly estimated, the origi- 
nal position of the sources in the HRD is fully recovered, 
rather, this indicates that ages are somewhat degener- 
ate in the sampled parameter range. In fact, as evident 
from Figure O upper panel, most of the stars are below 
the turn-off, and for the lowest masses, sources are lo- 
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TABLE 1 

Photometric filters included in TA-DA 



Photometric system 



Filter 



2MASS 
BESSELL 
CFHT 
COUSINS 
GALEX 
HST/ACS HRC 

HST/ACS SBC 
HST/ACS WFC 

HST/NICMOS 
HST/STIS CCD 
HST/STIS FUV 
HST/STIS NUV 

HST/WFC3 IR 

HST/WFC3 UVIS 



HST/WFPC2 



ISAAC 
JOHNSON 
JWST/NIRCAM 



LANDOLT 
SDSS 

SPITZER/IRAC 
STROMGREN 
TYCHO 
WFI 

WISE 



J, H, Ks 

U, B, BW, V, R, I 
I, Y, J, H, Ks 
R, I 

FUV, NUV 

F220W, F250W, F330W, F344N, F435W, F475W, F502N, F550M, F555W, F606W, 
F625W, F658N, F660N, F775W, F814W, F850LP, F892N 
F115LP, F122M, F125LP, F140LP, F150LP, F165LP 

F435W, F475W, F502N, F550M, F555W, F606W, F625W, F658N, F660N, F775W, 
F814W, F850LP 

FllOW, F160W, F165M, F187W, F190N, F205W, F207M, F222M 

50CC, F28X50LP 

25M, F25LY, F25Q, F25SRF2 

25M, F25CIII, F25CN182, F25CN270, F25MGII, F25Q, F25SRF2 

F098M, F105W, FllOW, F125W, F126N, F127M, F128N, F130N, F132N, F139M, 

F140W, F153M, F160W, F164N, F167N 

F200LP, F218W, F225W, F275W, F280N, F300X, F336W, F343N, F350LP, F373N, 
F390M, F390W, F395N, F410M, F438W, F467M, F469N, F475W, F475X, F487N, 
F502N, F547M, F555W, F600LP, F606W, F621M, F625W, F631N, F645N, F656N, 
F657N, F658N, F665N, F673N, F680N, F689M, F763M, F775W, F814W, F845M, 
F850LP, F953N, FQ232N, FQ243N, FQ378N, FQ387N, FQ422M, FQ436N, FQ437N, 
FQ492N, FQ508N, FQ575N, FQ619N, FQ634N, FQ672N, FQ674N, FQ727N, FQ750N, 
FQ889N, FQ906N, FQ924N, FQ937N 

F170W, F255W, F300W, F336W, F336W-no-leak, F336W-leak, F380W, F439W, 
F450W, F467M,F502N, F547M, F555W, F569W, F606W, F631N, F656N, F673N, 
F675W, F702W, F785LP, F791W, F814W, F850LP 
F1215, F1710, F2090, F3280, FLBB 
U, B, V 

F070W, F090W, F115W, F140M, F150W, F150W2, F162M, F164N, F182M, F187N, 

F200W, F210M, F212N, F225N, F250M, F277W, F300M, F322W2, F323N, F335M, 

F356W, F360M, F405N, F410M, F418N, F430M, F444W, F460M, F466N, F470N, 

F480M 

U, B2, S3, V 

u, g, r, i, z 

CHI, CH2, CH3, CH4 
u, Uoid, b, V, y 
B, V 

MB571, MB620, MB753, MB770, MB851, MB852, MB853, MB870, U, (7841, B, B842, 
V, R, I, 7879, Ha 
Wl, W2, W3, W4 
White throughput 



cated on the same main sequence regardless their age. 
Figure [7] confirms this: the error in the estimated ages 
(the difference between the best-fit age and the original 
age) is largest for the lowest masses in the sample, and 
progressively decreases towards the high-mass end. 

6. HOW TO OBTAIN TA-DA 

TA-DA can be retrieved from the following URL: 
http : //www. rssd. esa. int/ staff /ndario/TADA/. 

From this page the installation package, the manual, as 



well as additional grids of models can be downloaded. 



TA-DA is developed thanks to funding from the Na- 
tional Aeronautics and Space Administration (NASA) 
through award grant number NNX07AT37G. The au- 
thors thank Maddalena Reggiani, Abhijit Rajan, Marc 
Postman, and Jesiis Maiz Apellaniz for their help in the 
development of this software, and the team of scientists 
that proposed the project. 



APPENDIX 

EVOLUTIONARY MODELS, SYNTHETIC SPECTRA AND FILTER THROUGHPUTS AVAILABLE 



TA-DA includes the followin g evolutionary models: for PMS stars the iPalla &: Stahleri (199£), the 
iD'Antona fc Mazzitellil (|l998l ). the iBaraffe eFalT (fTOOl . the iSiess et all (|20Q0l) and the iTognelh'eraLl |20lT[ ). The 
first three of t h ese ar e available only for solar metallicity, whereas the last two include different values of [M/H]. The 
iTognelli et al.l ()2011[ ) models, furthermore, include isochrones and tracks for several assumed values of the mixing 
len gth parameter a, as well as different initial He and Deuterium abundances. For evolved population we integrated 
the lMarigo et al.l (jlOOS) isochrones. 

For what concerns the atmos phere models, TA -DA integrates the following grids: the [Kurucz (1993) models, the 
phoenix/NextGen models from iHauschildt et al.l (|1999D . the AMES-MT models from lAUard et al.. (■2000 ) (comple- 
ment ed with the NextG en grid for Tes > 5000 K), the AMES-Settl from |Allard_eLalJ (|2001|) , the BT-Settl models 
from lAllard et all (|2011f) . These grids cover different ranges of Toff, surface gravity logg, and metallicity. We refer the 
reader to the manual of TA-DA for these specifications. 
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Fig. 8. — A portion of the evolutionary tracks from'Si ess et al.1 l|2000l) . for the PMS-MS turn-on at intermediate masses. The left panel 
shows two isochrones, (1.6 and 2 Myr), obtained by connecting the points with the considered ages on each original track. The right panel 
shows the same isochrones obtained from the improved interpolation technique described in the text. 



TA-DA integrates a large number of band profiles for most of the commonly used photometry systems. These are 
listed in Table [1] we highlight that, although TA-DA is natively oriented to the optical and NIR wavelength range 
(since it addresses mostly problems linked to stellar photospheric fluxes), filters spanning from the UV to the mid- 
IR have been included. In the latter case (see, e.g, the Spitzer/IRAC and the WISE filters), the fact that TA-DA 
natively allows one to consider flux units (Jansky or erg s~^cm~^A~^) - commonly used at long wavelengths - besides 
magnitudes, simplifies the use of this tool. 

We stress that additional interior models, synthetic spectra, and filter throughputs can be added to the tool by the 
user with no need to modify the code. This possibility opens to future expansions of the capabilities of TA-DA. 

A NOTE ON THE INTERPOLATION OF PMS ISOCHRONES IN THE HRD 

It is not uncommon, especially for PMS models, to be computed with a sparse sampling of the stellar masses, and 
this can leads to a "segmented" appearance of the isochrones. This issue may become somewhat more serious in some 
ranges of the mass-age space, where the stellar evolution is very fast in the HRD and the rapidity of the evolution 
of Tcff and L changes quickly for contiguous mass tracks. This is the case of the PMS-MS turn-on, as illustrated in 
Figure m left panel. In this example we are showing the PMS tracks from lSiess et al.l ()2000D . as originally computed 
by these authors, and 2 isochrones for 1.6 Myr and 2 Myr, obtained in the standard way by connecting with straight 
segments the points, on all the tracks, corresponding to a given age. Between ^ 2.5 and 3.5 M0 the two isochrones 
cross, as it should not happen, and the segment connecting the 3.0 and 3.5 M© of the youngest isochrone passes in 
the region of the HRD between the 2.7 and 3.0 M0 tracks. 

Obviously, the ideal solution to this problem would be to have evolutionary models originally computed with a finer 
sampling of masses; when this is not possible one can only interpolate. The correct way to perform this is to intepolate, 
from neighboring mass tracks, points corresponding to identical evolutionary stages. This is the common practice, e.g., 
to derive post-MS isochrones, where the timescales of different evolutionary stages are very different between them 
and highly dependen t on the mass. For PMS stars, this method of interpolation is used, e.g., in the Pisa models from 
iTogneUi et all ((Ml . 

Unfortunately, this is not always possible. For example the iSiess et all (|2OOO0 PMS tracks are not provided for 
identical evolutionary points on different tracks. To improve on this, and compensate the issue shown in Figure |8l 
left panel, we propose an alternative method to perform this interpolation, maximizing the information present in the 
models. Our method is based on noting that, besides the 4 parameters listed in the evolutionary tracks (mass, age, 
logTeff, logL, contiguous mass tracks tend to have very similar shapes and sizes in the HRD, despite showing abrupt 
changes in the rapidity of evolution in the HRD. This is guaranteed by the fact that stars of similar mass undergo the 
same physical stages during their evolution, despite the fact that these may occur at very different times. Thus, we 
proceed at follows. First we interpolate new (interposed) mass tracks based in the shapes and lengths of the available 
ones. This can be done by assuming a metric in the HRD (e.g., a length I oc (10 logT)^ -f logL^), divide each pair 
of contiguous mass tracks into a fixed, large, number (n>1000) of segments of identical length, and linearly interpolate 
all the 4 parameters listed above between corresponding pairs of points on both tracks. Qualitatively, our approach 
may be thought as a "morphing" technique, which progressively "deforms" one track into the next one based on their 
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shape, size and position in the HRD. 

The result is shown in Figure [H right panel. We plot 20 intermediate mass tracks interpolated between the original 
3 and 3.5 M0 ones, computed as described. These now nicely follow the shape of the two original ones. We also 
show the same 2 isochrones as in the left panel, now obtained connecting points of same age on all the intermediate 
interpolated tracks. These isochrones now never cross, and also appear much s moother in the HRD. 

Some of the PMS models included in TA -DA (specifically, the Sie ss et al.l 120001 iBaraffe et all 119981 : 
iD'Antona fc Mazzitelli|[l998l : iPalla fc Stahleiill999f ) have been pre-interpolatcd using the approach we described, cre- 
ating very dense grids of masses and ages, generally including > 100,000 points per model. 
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